Homodyne measurement of average photon number 
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We describe a new scheme for the measurement of mean photon flux at an arbitrary optical 
sideband frequency using homodyne detection. Experimental implementation of the technique re- 
quires an AOM in addition to the homodyne detector, and does not require phase locking. The 
technique exhibits polarisation, frequency and spatial mode selectivity, as well as much improved 
speed, resolution and dynamic range when compared to linear photodetectors and avalanche photo 
diodes (APDs), with potential application to quantum state tomography and information encoding 
using an optical frequency basis. Experimental data also directly confirms the Quantum Mechanical 
description of vacuum noise. 

PACS numbers: 42.50.-p, 42.50.Xa, 42.50.Ar 



I. INTRODUCTION 

The most fundamental optical measurement is the in- 
tensity of a specific mode. Quantum optics tells us that 
the intensity is quantized |l|, appearing as a discrete 
photon flux. For a sufficiently small flux the individ- 
ual photons can be resolved and counted using a num- 
ber of standard techniques such as avalanche photodi- 
odes (APDs), photomultipliers and bolometers [2]. The 
optimum approach for a given application is dependent 
upon the wavelength and expected photon flux. An im- 
portant point is that at optical frequencies the thermal 
background is at zero temperature to an excellent ap- 
proximation. Thus "dark counts" , i.e. photon events in 
the absence of illumination, are technical in nature. 

It is also possible to probe an optical mode by first 
mixing it with a strong reference field, a local oscilla- 
tor, before detecting its intensity. This is called homo- 
dyne detection. In this case detection is generally by 
linear positive-intrinsic-negative (PIN) photodiode detec- 
tors which cannot resolve individual photons. Instead, 
the continuous spectrum of the signal's quadrature am- 
plitudes are measured. Quantum optics tells us that the 
quadrature amplitudes exhibit zero temperature fluctua- 
tions, thus a fundamental noise floor is observed for ho- 
modyne detection in the absence of signal illumination. 

Although all this is well known and the performance of 
photon counters and homodyne detectors has been tested 
individually for many systems 3] , to our knowledge a di- 
rect comparison of the photon counting and homodyne 
signal of the same field has not been made. This is per- 
haps why some authors claim that the homodyne statis- 
tics obtained in some continuous variable quantum infor- 
mation experiments are inconclusive |^. In this paper 
we make such a direct comparison. 

The paper is laid out in the following way. In section II 
we detail the theoretical relationship between continuous 
variable homodyne detection and the direct measurement 
of average photon flux densities. The experimental ap- 
proach used for homodyne detection is described in Sec- 



tion III and the results are compared with the measure- 
ments using a single photon detector (SPDM) in Section 
IV. Finally, a discussion of the significance of the results 
is given in Section V. 

As well as the fundamental interest of our results the 
demonstrated technique is shown to be easy to add to ex- 
isting experiments, providing fast measurements over a 
wide dynamic range. This potentially provides for signifi- 
cant reductions in the time taken to construct the density 
matricies for qudit systems [f| and other applications in 
quantum tomography. 



II. THEORY 

The annihilation and creation operators of an arbitrary 
field may be expressed in the Heisenberg picture as a sum 
of the amplitude X^ and phase X~^ quadrature opera- 
tors, where the annihilation A and creation A^ operators 
are given by pj 



A = 

it = 



x++ix- 



The quadrature operators represent continuous vari- 
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FIG. 1: Generic measurement concept, as applied to deter- 
mine n at the output of an unknown process. The direction 
of optical frequency translation as performed by the AOM is 
unimportant, as discussed in the text. 
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able observables that are able to be measured via homo- 
dyne detection techniques. 

If the field operator is expressed as a linear sum of 
steady state coherent amplitude A and quantum mechan- 
ical fluctuation terms 8A(t) such that A(t) = A + SA(t) 
and A(ty — A + SA(ty, and use is made of the Fourier 
transform relationship A(t) — > 8A(uo) and A(ty — > 
6A(— uj)^ , then the power spectrum of the fluctuations 
of an arbitrary input field can be written jy| : 

V ± {lo) = (5A(u)*5A(u) + 5A{-u)^5A{-u) 
±8A(-u)6A(u) ± 5A(-u)UA(u)^) 
+1 (1) 

where the spectral variance V^(u>) is given by 
(|5X ± (w)| 2 ). The first term of Eqn. is the average 
number of photons n + in the input mode at frequency 
the second represents n_ at — ui and the third and 
fourth terms represent correlations between the +uj and 
— u> sidebands. The final unit term represents the vac- 
uum fluctuations. In general, this may be expressed as 
0: 

where n(u>) = (n + + nJ)j1 represents the average num- 
ber of photons in both the positive and negative side- 
bands for a continuous variable measurement of the vari- 
ances at frequency to. Note that this relationship also 
remains valid for any pair of orthogonal measurements of 
the optical phase space, i.e. V e and V e+ % . In this case, 
V 9 {uo) = {\5X 9 {u)\ 2 ) where 6X e = (Sae w + Sa^- 10 ), 
representing an arbitrary quadrature operator at angle 9 
with respect to a specified phase reference. 

As a first example consider a coherent state \a). Co- 
herent states can be defined as eigenstates of the anni- 
hilation operator, i.e. A\a) — a\a). To a good approxi- 
mation such states can be produced by a well stabilized 
laser. We have V + = 4i?e(a) 2 + l and V~ = 4Im(a) 2 +l. 
Using Eqn.^Jwe find n = Re(a) 2 + Im(a) 2 = \a\ 2 which 
is as expected. 

As another example, consider the homodyne measure- 
ment of a squeezed state \a,r). Squeezed states can be 
defined as eigenstates of the operator B = y/G—WG — 1. 
Such a state possesses the following properties [lj: 

AX+ = e- r , AX- = e r 
(N) = \a\ 2 +sinh 2 r 

where the squeeze factor r = — ln(VG — y/G~—T). As- 
suming a similarly squeezed vacuum state (i.e. a = 0) 
with variances V + and V~ measured at frequency u>, 
Eqn. [3 yields: 



|e-'-| 2 + |e'-| 2 -2 

n (lj) — 

V ; 4 

= |0| 2 + sinh 2 r 

again in agreement with the expected result. The inter- 
esting thing in this second example is that the fluctua- 
tions in one quadrature can fall below the level associated 
with no illumination. This example shows that the vac- 
uum noise floor cannot be attributed to technical noise 
of the detector. 

The relationship of Eqn. [5] may also be exploited to 
measure the average photon number in a field of interest 
using the scheme illustrated in Fig. ^ In order to avoid 
low frequency technical noise, an acousto-optic modula- 
tor (AOM)y| introduces a radio frequency offset uo be- 
tween the mode of interest and a homodyne detector lo- 
cal oscillator. As shown, the mode of interest is upshifted 
relative to the local oscillator, and the vacuum mode in- 
troduced at frequency —uj relative to the local oscillator. 
From Eqn. □ (<5A(-w)t SA(-oj)) = (6A(-lj)SA(lu)) = 
(SA(— u>y5A(uiy) — and local oscillator phase depen- 
dance disappears. Consequently, V + = V~ and Eqn. |3 
yields: 

n+(uj) = (5A(u>y5A(u>)) = V{u) - 1 (3) 

where fi+(u>) is the average number of photons in the 
+lo sideband. A similar result is obtained if the +uj and 
— uj sidebands are interchanged, i.e. upshifting the local 
oscillator with respect to the field of interest. 

III. EXPERIMENT 

The experimental setup illustrated in Fig. [21 was an 
implementation of the topology shown in Fig. ^ The ex- 
periment was driven by a Photonetics Tunics- 1550 tun- 
able diode laser operated at 1540 nm, the wavelength 
chosen to permit stable mode-hop-free operation, as well 
as to take advantage of commercial telecommunications 
fiber products. The output beam was divided asymmet- 
rically, the majority of the power used as a coherent local 
oscillator for the homodyne detector and the remainder 
upshifted and attenuated to become the field of inter- 
est. The field of interest was then measured by three 
independent detection systems; a homodyne detector, an 
id-Quantique (id-200) InGaAs/InP SPDM and a NIST- 
traceable (Newport 3227) power meter. All interconnec- 
tions were made with single mode glass fiber (SMF) in 
the interests of reproducibility of measurement and for 
spatial mode filtering. 

To permit valid comparisons to be made between the 
three measurement techniques, all measurements were re- 
ferred back to Point (A) of Fig. [3 accounting for the 
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FIG. 2: Schematic representation of the system employed, 
modematching optics omitted for clarity. Points A, B\ and 
B2 represent the states to be measured. 



losses unique to each optical path in the analysis of re- 
sults. As physical disturbance of the SMF altered the 
polarisation of the propagating photons via the stress- 
optic effect 9], a series of data was recorded using a 
50/50 fiber coupler to permit simultaneous SPDM and 
homodyne measurements to be made whereby the signal 
path waveplates had been previously adjusted for maxi- 
mum homodyne efficiency. This also made long time con- 
stant drifts in the laser's output power common to both 
measurement systems. However, to support the claim 
that both measurement schemes were truly independent, 
a small set of results were collected by interchanging the 
connections of one port of the coupler between the SPDM 
and the homodyne detector. 

Experimentally, measurement of the mean photon 
number is determined by integrating the mean photon 
flux $ over time period r. The field of interest at point 
(A) may be approximated as a weak coherent state \(x)ai 
denned as 



\a) = e 2 



£ 

71=0 



(4) 



where the mean photon flux is related by $ = \a\ 2 = 
PX/hc where P is the optical power, A is the optical 
wavelength, h is Planck's constant and c is the speed of 
light. 

\a) a is further attenuated in a calibrated fashion to 
yield |a) b 12 with o>b 1 2 varied by adjusting the RF drive 
power to the AOM to yield values for 2 in the 
range ~ 10 3 — ► 10 5 photons/sec. The second order 
Bragg-diffracted optical output was used to achieve 
a greater attenuation than was otherwise possible with 
the first order output, and to avoid radiative coupling 



between the RF signal generator and the homodyne de- 
tector. Hence, the AOM was operated at its frequency 
of peak diffraction efficiency of 80 MHz, and homodyne 
detection performed at a center frequency of 160 MHz. 

The homodyne detector was implemented with 
matched InGaAs PIN photodetectors and an Agilent 
E4407B spectrum analyser (SA). Free space optical com- 
ponents were used to permit independent control of in- 
put polarisation and relative phase, the latter via a PZT 
mounted mirror. 

As the standard deviation a n of the Poissonian distri- 
bution described in Eqn.^is related to the mean photon 
count n by ofj = h, the variance of the SPDM measure- 
ments is reduced by maximising r. Integration times of 
5 minutes with an APD gating period tq of 100 ns at 
10 /is intervals were used to satisfy the practical con- 
straints imposed by long time constant laser power drifts 
and afterpulsing probability [lljj]. The dynamic range of 
the SPDM measurements were restricted to only three 
orders of magnitude, the upper limit imposed by satu- 
ration of the APD gated time bins and the non-linearity 
associated with the increasing probability of > 2 photons 
arriving within each gate window. The lower limit arises 
as a consequence of the finite dark count probability and 
low observed SPDM quantum efficiency r\s of 11%. 

Homodyne measurements were performed with a nom- 
inal SA resolution bandwidth (RBW) of 30 Hz and mea- 
sured RBW of 33.18 Hz. This corresponds to an effective 
integration time of t/RBW = 30.1 ms, data collected 
for a 30 s period. The overall efficiency of the homo- 
dyne detector rj is a function of the quantum efficiency 
rjdet of the matched PIN photodetectors and the inter- 
ferometric fringe visibility VIS, where rjn — VIS 2 and 
rj = rjHijdet- Accounting for the transition from fiber to 
free space, the maximum experimental fringe visibility 
attained was 93%, the corrected mean photon number 
observed by the homodyne detector thus given by Eqn. 
El scaled by RBW/rj. 



IV. RESULTS 

Fig.Hillustrates the SPDM inferred mean photon flux, 
showing a useful linear relationship over approximately 3 
orders of magnitude, with the measurement uncertainty 
increasing rapidly at low light levels. Following the cor- 
rection of technical noise as detailed in the following sec- 
tion, data points with a negative mean have been deleted 
and error bars truncated at $ = 0. 

Note that for both plots, the incident photon flux was 
determined as per Fig. [2] via macroscopic power mea- 
surements and the solid line represents the ideal 1:1 rela- 
tionship between incident and measured values. The er- 
ror bars for both axes are placed one standard deviation 
from the mean. Data points obtained via the coupler are 
marked with a point, the circled points (seen at incident 
fluxes of 1.6 x 10 5 , 8.0 x 10 6 and 9.8 x 10 8 photons/sec.) 
represent data obtained by interchanging connections as 
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Incident photon flux (photons/sec) 



FIG. 3: Average photon number at point (A) of Fig. [5] as 
determined by SPDM measurements. 




Incident photon flux (photons/sec) 



FIG. 4: Average photon number at point (A) of Fig. as 
determined by homodyne detection techniques. 



detailed in the proceeding section. 

The homodyne detector results are shown in Fig. 01 the 
point /circled data values determined via Eqn. [3] exhibit- 
ing a comparatively superior dynamic range in excess of 
4 orders of magnitude. If the subtraction of the vacuum 
noise contribution is neglected as required by stochastic 
electrodynamics then Eqn.^lreduces to n(uS) — V(uj), the 
consequence of this alteration illustrated by the diamond 
data points. Is is observed that whilst the magnitude of 
the uncertainties are reduced, the results clearly disagree 
with both the SPDM data in Fig. and the desired 1:1 
relationship. Although the two sets of results converge at 
higher flux levels where the magnitude of vacuum noise 
ceases to be significant compared with the measured vari- 
ance, low flux levels are only calculated correctly by the 
quantum noise model. 



As suggested in Eqn. [2 no sensitivity to relative optical 
phase was observed in the homodyne variance as the PZT 
illustrated in Fig. |21 was swept over a number of optical 
wavelengths. 



V. DISCUSSION OF RESULTS 

Other than the obvious differences in dynamic range, 
the demonstrated homodyne measurement technique dif- 
fers from the SPDM photon-counting approach in a num- 
ber of other aspects. 

Most significantly, the homodyne detector is a con- 
tinuous variable measurement system, and thus is un- 
suitable for conditioning purposes, such as required by 
non-deterministic qubit operations for implementation of 
linear optical quantum computation |12fl . 

Homodyne detection is also inherently sensitive to the 
spatial mode, polarisation and frequency of the input 
state to be measured. The RBW of the SA dictates the 
frequency selectivity of the technique, with the maximum 
detection bandwidth determined by the RF components, 
i.e. the linear photodetectors, subtractor and SA. As de- 
tection is reliant on optical inference at the beamsplitter, 
the spatial and polarisation mode of the LO and signal 
beams must be identical. This provides a high degree 
of immunity to the detection of ambient light as such 
photons are incoherent with either the signal or LO. 

Conversely, the SPDM is insensitive to the input spa- 
tial mode, and exhibits frequency sensitivity only in re- 
sponse to the bandwidth of the input optics and the spec- 
tral responsivity of the semiconductor APD. Experimen- 
tally, this gave rise to two sources of systematic error; the 
coupling of ambient room light to fiber cladding modes 
which could be controlled by shielding the SPDM SMF 
input, and that due to the AOM isotropically scattering 
a portion of the (unshifted) input beam and this light 
being modematched into the SMF. Reducing the AOM 
optical input power by changing the preceding beamsplit- 
ter ratio reduced the count rates arising from this effect 
to a level comparable to the APD dark count rate. 

Technical noise adversely affects both measurement 
schemes. In the SPDM, technical noise manifests itself as 
undesired count events and missed detection events, these 
effects becoming dominant at low flux levels. This is pri- 
marily a consequence of low quantum efficiency and dark 
counts, both of these a function of the APD semiconduc- 
tor material. The time series for the SPDM shown in Fig. 
[5] illustrates typical data when the SPDM is illuminated 
(upper trace) and in the absence of optical input (lower 
trace). The probability density function of the "bright" 
data is given by Eqn. 0] and the distribution for the dark 
counts is also Poissonian in nature 01 . As both distri- 
butions are of the same family, the mean of the technical 
noise may simply be subtracted away [l4| to yield the 
mean number of detection events, provided that the two 
distributions may be resolved. 

In the homodyne scheme, technical noise is comprised 
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FIG. 5: SPDM time series and statistical distributions for $a ~2x 10 7 photons/sec. 
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FIG. 6: Homodyne detector time series and statistical distributions for <3>a 
distributions associated with the technical noise sources are illustrated. 



2 x 10 7 photons/sec. For clarity, only the 



mainly of electronic noise, the primary sources being the 
SA noise floor and dark noise from the detectors, i.e. 
noise generated in the absence of optical input. The 
time series data and associated log-Rician 0, 0] dis- 
tributions for these sources are shown in Fig. |SJ along 
with the detected quantum noise limit (QNL) and de- 
sired signal. Although, as above, the contribution of dark 
noise may be subtracted away from the measured QNL, 
the noise serves to decrease the precision of the measure- 
ment. This deleterious effect may thus be minimised by 
increasing the separation between the observed QNL and 
the technical noise. As the detector photocurrent linearly 
scales with the optical input intensity within limits im- 
posed by saturation Q , the increase in signal to technical 
noise ratio may be readily achieved by raising the local 



oscillator power. 



The statistical uncertainty in each measurement oc 
TV -1 where N is the number of independent sample 
points considered |16|| . As dominates the final un- 
certainty in the measurement precision may be im- 
proved by increasing N within practical time constraints. 
As demonstrated, the homodyne detector permits much 
larger values of N to be attained whilst the total inte- 
gration time Nt = 30 s remains small compared to 250 
s for the SPDM. 
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VI. CONCLUSION 

We have experimentally verified the proposed homo- 
dyne measurement technique and demonstrated a supe- 
rior sensitivity, dynamic range and measurement speed 
than is otherwise possible via linear PIN photodetectors 
or APD based SPDMs. The homodyne detection scheme 
also offers intrinsic polarisation, spatial and frequency 
selectivity, and does not require elaborate phase locking. 

The results of simultaneous quadrature and APD mea- 
surements also clearly differentiate between semi-classical 
and fully quantum models of optics. In particular, semi- 
classical models such as stochastic electrodynamics 01 



are unable to account for the observed quadrature vari- 
ance at the QNL in the absence of photon detection 
events of the SPDM Q. 
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